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ABSTRACT 

Aims. We present an update of the 3C 273 's database hosted by the ISDC, completed with data from radio to gamma-ray observations 
over the last 10 years. We use this large data set to study the multiwavelength properties of this quasar, especially focussing on its 
variability behaviour. 

Methods. We study the amplitude of the variations and the maximum variability time scales across the broad-band spectrum and 
correlate the light curves in dilferent bands, specifically with the X-rays, to search for possible connections between the emission at 
dilferent energies. 

Results. 3C 273 shows variability at all frequencies, with amplitudes and time scales strongly depending on the energy and being 
the signatures of the different emission mechanisms. The variability properties of the X-ray band imply the presence of either two 
separate components (possibly a Seyfert-like and a blazar-like) or at least two parameters with distinct timing properties to account for 
the X-ray emission below and above ~20 keV. The dominant hard X-ray emission is most probably not due to electrons accelerated 
by the shock waves in the jet as their variability does not correlate with the flaring millimeter emission, but seems to be associated 
to long-timescale variations in the optical. This optical component is consistent with being optically thin synchrotron radiation from 
the base of the jet and the hard X-rays would be produced through inverse Compton processes (SSC and/or EC) by the same electron 
population. We show evidence that this synchrotron component extends from the optical to the near-infrared domain, where it is 
blended by emission of heated dust that we find to be located within about 1 light-year from the ultraviolet source. 

Key words. Astronomical data bases: miscellaneous - Galaxies: active - quasars: individual: 3C 273 - X-rays: galaxies - Infrared: 
galaxies - Ultraviolet: galaxies 
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1. Introduction of an e xternal X-ray source (ICollin-Souffrinlll991h or of a hot 

S/^ corona (lHaardtetalJll994 . and seems to be due to the con- 

Active galactic nuclei (AGN) are known to emit a broad-band tribu tion of two compon ents with different variability proper- 

spectrum that extends from radio to gamma-rays. Our under- jjes dPaltani etalJI 1998^ 1. In the X-rays, it has been suggested 

standing of their multiwavelength emission relies on the iden- that the soft-excess could be due to t hermal Comptoni sation of 

tification of all components that contribute to the radiation. cool-disc photons in a warm corona (iPage et alJl20"04h . Inverse 

3C 273, a bright and nearby (z = 0. 158) radio loud quasar, is a Compton processes of a thermal plasma in the disc or in a corona 

well known example of an AGN whose emission has been ex- ^nd of a non-thermal plasma associated to t he jet are believed to 

tensiv ely studied in the past at all frequencies ( CourvoisieretalJ generate the X-ray to gam ma-ray emission (iKataoka et alJl200l 

119871 1 1990) . It shows most of the properties characteristic of iGrandi & Palumbol2004l) . 

blazars, Uke strong radio emission, a jet with apparent superlu- ^he 3C 273's Database hosted by the ISDCQ is one of the 

mmal motion, large flux variations and, occasionally, polarisa- ^^^^ complete multiwavelength database s currently availabl e 

tion of the optical emission (for a review, see|Courvoisier 1998). ^GN. When it was first published (iTtirler et al.lll99"9bl) . 

Synchrotron flares from a relativistic jet dominate the radio-to- contained 70 light curves covering 16 orders of magnitude in 

milhmeter energy output and extend up to the infrared and op- ^^^^^^ ^^^^gy 30 y^^^ observations. We present here 

tical domams dRobson et al. || 1993 |; i Tiirler et alj | 2000|) , whereas ^„ ^p^^^g ^f the database with recent data and special emphasis 

thermal emission from dust is at least m part responsible for the j^e X-ray emission. Using these data, we analyse the multi- 

mfrai-ed quiescent continuum (RobsonetalJll983; TurleretalJ wavelength variability properties of 3C 273, studying both the 

12006). The particulai-ly bright excess in the optical-UV band has amplitudes and the time scales of the variations as a function of 

been interpreted as a signature of the accretion disc iShieldsl j^e frequency and ti-ying to understand the different variability 

[19781 and many subsequent studies), maybe in the presence behaviours in the framework of already existing emission mod- 
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Fig. 1. Time coverage of all observations of the database as a function of the frequency. The colour-scale indicates the intensity of 
the observed flux in vFy representation and spans from a minimum intensity of 6 x 10"''* to a maximum of 1 .6 x 10"' erg cm"^ s"' . 



els. We cross-correlate the fluxes in different energy ranges to in- 
vestigate the connections between the different emission compo- 
nents, especially concentrating on the origin of the X-ray emis- 
sion. 



2. Data collection 

We included all the 3C 273 public data of the last 1 years in the 
databas e (Table [TJ, following the same criteria as iTiirler et aU 
(Il999bl hereafter Paper I). Fig.[T]shows the time coverage of the 
data as a function of frequency. 



2. 1. Radio and millimeter data 

In the radio, we added the recently published data at 22 and 
37 GHz from the Metsahovi Radio Observatory (Teras ranta et al.l 
|2004 l2005h . complemented by unpublished data from this ob- 
servatory and from the University of Michigan Radio Astronomy 
Observatory (UMRAO) at 4.8 GHz, 8.0GHz and 14.5 GHz. 
We also added the measur ements from the lOOn i antenna at 
Efelsberg published by Reich et al.l (Il998l) and iTiirler et alj 
(|2006|) . Finally, we completed some historic Ught curves with 
early data not published in tabular format by extracting them 
from figure s. This was done for the earliest Haystack data at 
15.5 GHz inlAllen et alj(|l968|), the UMRAO data at 8.0 GHz in 



lAUer & Haddock (,1967 ) and some Crimean time-averaged data 



EfanovetalJ (119771) . 



at 22 GHz derived from 

At millimeter and submillimeter (mm/sub-mm) wavelengths 
we added the observations published by Robson et al. (200l|) 
from the SCUBA instrument at 850/im on the James Clerk 
Maxwell Telescope (JCMT), Hawaii and more recent unpub- 



lished observations by this telescope. Other mostly unpub- 
lished observations are from the Swedish-ESO Submillimeter 
Telescope (SEST) in La Silla, Chile, from the Submillimeter 
ArrajQ (SMA) on Mauna Kea, Hawaii (for information on the 
quasar monitoring program at the SMA see Gurwell et al. 2007j), 
the Australia Telescope Compact Array (ATCA), and the Owens 
Valley Radio Observatory (OVRO), California. 



2.2. IR data 

An important addition t o the infrared dataset are the lon g-term 
light curves published bv'Neugebaue r & Matthewsid 19991) in the 
J, H, K, L and N bands from the Palomar Hale 5 m telescope. 
Unfortunately, these data could not be obtained in tabular form 
and have therefore been extracted from the figure and converted 
into flux densities usin g the zero-magnitude fluxes given by 
'Neugeba uer et al.l (11979 ). Other interesting addi tions are the data 
obtained during synchrotron outbur sts by McHardv et al.l d 19991 
120071) . the data of April 2002 fromlJorstad et al.l (|2004 and the 
data of June 2004 from ITiirler et all (l2006h . In the far-infrared, 
we added the obse rvations by the In frared Space Observatory 
(ISO) pubUshed bv lHaas et al.l(l2003h . 

The European Southern Observatory (ESO) magnitudes are 
now conv erted to flux using the zero-magnitude fluxes pub- 
lished by iBouchet et al.l (Il99l]) . To have the most coherent 
dataset with the observation of the United Kingdom Infrared 
Telescope (UK I RT), w e now also rescaled the ESO data from 
iLitchfield et al] (1 19941) based on the ratio of ESO and UKIRT 
zero-magnitude fluxes. Other data than those mentioned above 



http://smal . sma. hawaii. edujtools. html 
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have all been converted from magnitude to flux with the UKIRT 
zero-magnitude fluxes that was used for all data in Paper I. 

2.3. Optical and UV data 

In the optical domain we included recent unpublished data in 
the Geneva photometric system using the Mercator telescope on 
La Palma, Canary Islands. We also added the observat i ons of 
the Wise Observatory, Israel published by iGiveon et alJ d 19991) 
and lKaspi et alj (l200Cft. Based on 17 quasi-simultaneous (within 
5 days) observations in the Geneva photometric system, we de- 
rived a scaling facto r of 1. 106, applied to the B-band flux den- 
sities of iKaspi et al .' (200(f). In addition, we included the I and 
R-band data from Jorstad et al. (2004) and the R-band data from 
lGhoshetal.l(l2000l) . 

The scaling of other data sets already inclu ded in Paper I did 
not change, except for the Burkhead system (lBurkheadlll980l 
and references therein) V-band magnitude correction that is now 
0.00 instead of -1-0.05. This difference comes from the correction 
of a quite obvious typo in the publication of Lvutvjj (119761) : B- 
V should be +0.20 instead of -t-0.02 in Hne 8 of Table 1. In the 
R and I bands, we no longer use the zero-magnitude fluxes of 
Paper I, but make a distinction, as far as possible, between data 
obtai ned in the standard Johnson system (R: 2941 Jy; I: 2635 Jy) 
or flie lCousinsI (119761) svstem (R: 3080 Jy; I: 2550 Jy). In partic- 
ular we used t he Co usins values for the R-band magnitudes from 
IGiveon et al.ld 19991) . 

We also included V-band measurements by the Optical 
Monitoring Camera (OMC) aboard the INTEGRAL satellite 
taken from the OMC Archival We averaged contiguous obser- 
vations and converted magnitudes to fluxes in the Geneva pho- 
tometric system as described in Paper I, without any additional 
correction (Vcorr = , see Paper I). 

Optical and UV data from XMM/OM observations have 
been taken from IChernvakova et all (l2007h . when available, or 
extracted following the procedure reported there. The 45 OM 
fluxes, covering from 7 to 12 different epochs, depending on the 
wavelength, have been included in the database after checking 
that there are no calibration problems and that the wavelengths 
of the 6 OM filters are sufficiently close to those of the light 
curves in the database. 

Among the few observations performed with the 
SwiftfUVOT monitor, only data collected with the V filter 
are included in our database, because, due to the brightness 
of 3C 273 (m = 11-13 depending on the filter), the detector 
is operating above its saturation limit with all filters except V. 
The V fluxes have been obtained analysing the data with the 
XRTDAS software v.2.5a distributed with the HEAsoft 6.1.1 
package and using the latest calibration files (2007-03-30) 
available at the time of the work. 

Finally we included in the database an additional light curve 
in the far-UV at a wavele ngth of 1050 A w ith data extracted from 
lAppenzefler etal] (Il998h and lKriss et al] (fl999.) . 



2.4. X-ray data 

During the last 10 years, several X-ray satellites have exten- 
sively observed 3C 273, considerably increasing the amount of 
data available above 0.5 keV (Fig. [T]i. In order to obtain a data 
set as uniform as possible, we fitted the available spectra in 
similar energy ranges using a power law model absorbed by a 



Galactic column density A^h = 1.8 x 10^" cm"^. Some simul- 
taneous or quasi-simultaneous observations (within about one 
day) of 3C 273 have been performed in the X-rays allowing 
the inter-calibration of the different instruments. Chandra ob- 
servations of 3C 273 were mainly conce ntrated on the X-ray jet 
(ISambruna et al.ll200ll; iJester et d]|2006l) and therefore have not 
been included in the database. 



2.4.^ . XMM-Newton/PH data 

Among the numerous observations of 3C 273 performed by 
XMM, 16 sets of EPIC/PN data collected from June 14, 2000 
to July 10, 2005 were publicly available at the time of this w ork. 
The EPIC spectra have b een taken from IStuhlinpeil (|2004|) and 
from IChernvakova et al.l (12007). Due to their higher signal-to- 
noise ratio, we choose PN spectra instead of MOSl/2 for the 
spectral analysis. The spectra are fitted separately in the ranges 
0.5-2 and 2-10 keV in order to obtain fluxes and spectral indices 
in the same energy bands as for BeppoSAX data. 



2.4.2. RXTE data 

Several campaigns have been performed with RXTE to monitor 
the X-ray emission of 3C 273. We have collected all the data 
available at the time of this work in the RXTE archive, which 
consist in 957 observations performed between February 9, 1996 
and December 11, 2005. We use the RXTE standard products to 
derive a flux and a spectral index associated with each obser- 
vation. The HEXTE data show a low signal to noise ratio and 
consequently a large scatter in the resulting parameters when fit- 
ted alone. Therefore we perform a PCA-HEXTE combined fit 
for each observation, resulting in a broad-band spectrum from 3 
to 100 keV. Our choice is supported by the absence of strong fea- 
tures, like Compton reflec tion, in the X-ray spectrum of 3C 273 
(iGrandi & Palumboll2004 . 

A slow gain drift of about 1 % over 2 yearfl determines small 
changes in the energy-channel relation of the PCA instrument. 
Consequently, the maximum energy of channel 2 moved from 
2.02 to 2.96 keV over the mission time. Since the inclusion of 
this channel results in very high reduced chi-square values xl^.^^ 
we exclude the data below 3 keV from all PCA spectral fitting. 

We choose XMM/PN as the reference instrument and renor- 
malise the RXTE fluxes to the PN ones at 5 and 10 keV, obtaining 
an inter-calibration factor of ~ 0.78. The light curves of the other 
X-ray instruments are then reseated either to the ZMM/PN one 
or to the RXTE one, depending on the time frame and on the 
energy band. The inter-calibration factors reported in Table |2] 
for BeppoSAXI'PYiS, INTEGRAL and SwiftlXRT are relative to 
RXTE, therefore they are multiplied by 0.78 to derive the final 
scaling factors (relative to XMM/PN) that we apply for the light 
curves in the database. 



2.4.3. Beppo-SAX data 

The LLCS, MFCS and PDS spectra of the 9 BeppoSAX observa- 
tions of 3C 273 have been retrieved from the ASDC archiv^ and 
separately fitted in the ranges 0.5-2 keV, 2-10 keV and 15-200 
keV, respectively. This allows us to have a separate representa- 
tion of the soft-excess region, the medium and the hard X-ray 
bands. The 0.5-10 keV fluxes are renormaUsed to the XMM/PN 



http://sdc. laeff. into, es/omc/ 
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Table 2. Inter-calibration constants for the X-ray light curves 



Energy [keV] 


RXTE 


Beppo5AX/LECS 


BeppoSAXIMECS, 


BeppoSAXrPm 


XMM/PN 


INTEGRAL 


5w//f/XRT 


0.5 


- 


1.67 ± 0.07 


- 


- 


1 


- 


4.87" 


1 




1.35 ± 0.02 






1 




4.87" 


2 






0.87 ± 0.02 




1 




4.87° 


5 


Q.ll ± 0.12 




1 




1 


1 


4.87* ± 1.27 


10 


0.78 ± 0.15 




1 




1 


1 




20 


0.78'- 






1 




1 




50 


0.78" 






1 




1 




100 


0.78" 






1 




1 




200 








1 




1 




Renormalised to 


XMM/PN 


XMM/PN 


XMM/PN 


RXTE 




RXTE 


RXTE 



" Fixed to the value found at 5 keV. 

* This large scaling factor is due to the pile-up affecting the XRT images of 3C 273 and it does not have to be considered as an inter-calibration 
constant between SwiftfXKI and RXTE. 
" Fixed to the value found at 10 keV. 



ones. The inter-calibration factor between the 20-200 keV fluxes 
from PDS and the RXTE ones is found to be consistent with 1 . 



2.4.4. INTEGRAL data 

From the beginning of the mission, INTEGRAL spent about 2 Ms 
on 3C 273, more than 1 Ms of which (January 2003-July 2005) 
are publicly available and included in the database. Data analy- 
sis is described in detail in Chernvakova et al.l (120071) . Different 
observations were summed up in order to increase the statistics 
of the data and resulted in 5 sets of IBIS/ISGRI and JEM-X data 
and 3 of SPI data. For each data set a combined fit with all the 
available high-energy INTEGRAL instruments is performed in 
the range 5-200 ke V. SPI (or ISGRI when SPI data are not avail- 
able) is assumed to be the most stable INTEGRAL instrument 
and the o ther inter-calibratio n factors are fixed to the values re- 
ported bv Chernv akova et alj (12007). The inter-calibration factor 
between INTEGRAL and RXTE fluxes is consistent with 1 . 



2.4.5. SwiftlXRl data 

Swift has observed 3C 273 several times, collecting data with all 
3 instruments on board. We have analysed and included in our 
database the 5 XRT observations performed in Photon Counting 
readout mode with exposure time greater than 3 ksec (between 
November 15 and December 19, 2005). The data were analysed 
with the XRTDAS software v.2.5a and the 2007-03-30 calibra- 
tion files. During all these observations the source shows a count 
rate higher than 1 count/s and the image is affected by photon 
pile-up. In order to correct for this effect, the spectra are ex- 
tracted in an annu lar region of 6 and 2 pixels inner and outer 
radii, respectively (IVaughan et al.ll2006l) . Therefore, only a frac- 
tion of the flux can be recovered, which is estimated thanks to 
quasi-simultaneous RXTE observations (Table|2]i. For December 
19, 2005 observation, we performed an exposure map coiTection 
on the data (using xrtexpomap) in order to account for the de- 
creased collecting area due to dead columns of pixels crossing 
the source image. The effect is found to be negligible, showing 
no changes in the fitting parameters compared to the uncorrected 
data. The fitting on all the XRT spectra is separately performed 
in the 0.5-2 and 2-8 keV energy bands. 



2.4.6. BATSE data 

In the hard X-ray range, the BATSE data from the Compton 
Gamma-Ray Observatory (CGRO) have recently been made 
available by the BATSE Earth occultation tearrO for four dif- 
ferent energy bands: 20^0, 40-70, 70-160 and 160^30 keV. 
In order to make them fit the database format, we extracted four 
different BATSE light curves for energies of 20, 50, 100 and 
200 keV assuming a fixed power law photon index of a - 1.7. 
No rescaling factor has been applied to the BATSE light curves. 
The variability analysis was performed on the 20-70 and 70^30 
keV light curves separately from the other hard X-ray data (see 
Sect.0. 

2.5. Gamma-ray data 

In the gamma-ray domain, we included the combined 3 MeV- 
to-lOGeV power-l aw fits to the CGRO/C OMPTEL and EGRET 
data pubhshed by ICollmar eTaP (l2000bh from 1996-1997 ob- 
servations, as well as single-energy-bin COMPTEL data from 
ICoIlmar et all (l2000a"F). Additio nal CGRO/EGRE T data come 
from spectral fits publis hed by iMukheriee et al.l (11997 ) and 
iNandikofliur et al.1 i2007l) . Finally, at 100 MeV, we included 
the EGRET fluxes at energies E > 100 MeV given by 
[Mukherjee et al. (1997) from observations in 1991-1994. We 
note that there are sometimes several entries in the database for 
the same observations analysed by different authors that can ei- 
ther add or reduce intrinsic source variability. 

3. Database format 

We took the opportunity of the update of 3C 273's database 
to reorganise slightly the data in different light curve files and 
to include additional information. The database, available at 
http://isdc.unige.ch/3c273/, is organised in 70 ASCII files, each 
of them containing a light curve at a certain wavelength. Each 
file is composed by 10 columns with one row per observation. 
Most of these columns (date in decimal year, frequency, wave- 
length, flux, flux error, reference) have maintained the same for- 
mat as in Paper I. The date in decimal year is still linearly re- 
lated to the Julian Date, but we now also give the date in UTC 
format as yyyy-mm-ddThh:mm instead of the previous format 
yyyymmdd . hh. For observations lasting several days, we express 

* http://f64. nsstc. nasa. gov/batse/occultation/ 
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the date as a range of the form yyyy-nun-dd/yyyy-imn-dd in 
UTC format and with the decimal year corresponding to the mid- 
dle of the observation. In some cases when the day of the month 
or even the month is not known we cut the UTC format to the 
last relevant number and assume the middle of the month (day 
15) or of the year (July 1st) to derive the decimal year 

To allow an easy selection of the data, we introduced a 
"Flag" column giving information on the quality and reliability 
of the data. It can have integer values from -3 to + 1 . The default 
value for good data points is zero and less useful or reliable data 
points have negative flag values. A flag of -1 indicates "useless" 
data in the sense that they are compatible with contemporaneous 
observations with smaller uncertainties, but their inclusion just 
adds noise to the light curves. A flag of -2 indicates "uncertain" 
data that are less reliable because of large uncertainties or due 
to an automatic processing without manual checking. Finally, 
a flag of -3 is assigned to "dubious" data points that are clearly 
outliers and incompatible with contemporaneous observations or 
with the overall variability trend. Finally, we assign a flag of +1 
to good infrared and optical data with a significant contribution 
by synchrotron flares lasting several days to weeks. 

The flagging of the data was usually done manually by care- 
ful inspection of the light curves. An automatic procedure was 
only used for the daily radio measurements from the Green 
Bank Interferometer (GBI). This was done for points that were 
outstanding by more than 3-cr from the 30-day smoothed light 
curve as we do not expect true variability on a much faster 
time-scale than a month at these low radio frequencies. For 
these data, as well as for the Earth occultation measurements 
by CGRO/BATSE in the hard X-rays, we both include the sin- 
gle, daily observations with a flag set to -1 and the more easily 
usable 30-day averages with a flag of zero. 

The meaning of the Flag value and additional information 
is included with each measurement in a new column called 
"Comments". For instance, we included there the scaling fac- 
tor used for the G BI data according to the values given by 
iTiirler et"an ( l2000l) . as well as information on the assumed hy- 
drogen column density A^h (fixed or free) for X-ray data from 
the literature. Finally, the information on the telescope or the 
antenna in the "Observatory" column was also extended in par- 
ticular by the addition of its diameter in meter 

In addition, at the same web page, for each light curve a fig- 
ure in GIF format is now provided, where the data points with 
different flag values are indicated by different colours. 
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Fig. 2. Examples of light curves from the 3C 273 database 
for the last 23 years of observations, at 5 GHz, 37 GHz, 
0.8 mm, in the K band, in the V band, at 5 keV and in 
the 20-70 keV range (this latter rebinned to 1 -month bins). 
The data during strong synchrotron flares (Flag =1) are in- 
dicated in grey (red in the electronic version) for the op- 
tical and IR data sets [for a higher resolution figure, see 
http://isdc. unige. ch/~soldi/image/9947F2Jir. ps\ . 



4. Variability analysis 

In order to use only the best-quality data for the variability anal- 
ysis, we consider only Flag > data and exclude the light curves 
with less than 10 data points. Moreover, due to inter-calibration 
problems, the radio light curves at v < 1.5 GHz and the I and R 
optical data sets are not included. 

Synchrotron flaring e mission is known to con tribute in the 
infrared-optical domain jCourvoisier et alj Il988l) . In order to 
separately study the variability behaviour of the underlying com- 
ponents, infrared and optical light curves excluding the peri- 
ods of known synchrotron flares (January-June 1972, December 
1982 - June 1983, Januai-y-May 1988, February-August 1990 
and January-March 1997, flagged with 1 in the database; see 
Fig.|2]i are also analysed ('wf data sets in Table[T]l. 

The 20-70 and 70-430 keV BATSE hght curves are used 
separately from the other X-ray data, as they have been rebinned 
to about one-month bins in order to have an acceptable signal- 
to-noise ratio for the variability analysis (Fig.|2]). 



The RXTE/PCA fluxes included in the database have been 
extracted from the spectra accumulated during observations 
from 0.2 to 17 ks long. In order to study variability on scales 
as short as possible we use the PCA light curves (in counts/s) 
provided by the RXTE archive in the 4-9 and 9-20 keV ranges, 
after rebinning them to 1000-second bins. 

In Table [T] the light curves are reported together with the 
time range between the first and last observation included in the 
database, the number of data points, the average frequency, the 
average flux with its standard deviation, the vFy intensity and 
the variability parameters that will be discussed in the following 
sections. Some examples of light curves from the database are 
shown in Fig.|2l The average spectral energy distribution (SED) 
including all the light curves of the database is shown in Fig. [3] 
Beside the synchrotron and Compton peaks, typically observed 
in the SED of blazars, 3C 273 shows a prominent excess in the 
optical-UV band, the so-called blue bump. 
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4.1. Fractional-variability amplitude 

4.1.1. Definition 

AGN are known to be very variable objects at all wavelengths 
and different methods can be used to quantify the amplitude of 
the variations. In order to estimate the total variability of each 
light curve in our database we use the fractiona l variability am- 
pHtude, Fvar, given by (e.g. lEdelson et al.|[2002h : 

where 5 ^ is the sample variance of the light curve, x is the av- 
erage flux and Seir - jjYji ^en , mean of the squared mea- 
surement uncertainties. Fvar is calculated for all the light curves 
in our sample (Table[T]) and the resulting spectrum, from radio to 
gamma-rays, is shown in Fig.|4] The Fym values obtained when 
excluding synchrotron flares from IR and optical light curves are 
also calculated, though in the optical range the flares are too few 
to have any significant influence on the value of Fyai - 

The uncertainty on Fyai due to the measurement-error fluc- 
tua tions have been estim ated through Monte Carlo simulations 
by Vaug han et"a!1 (l2003h (see Eq. (B2) there). Applied to our 
data sets, it results in percentage uncertainties smaller than 5% 
in more than 80% of the light curves. In order to give a more 
conservative estimate of the uncertainty we use a bootstrap pro- 
cedure resulting in error bars that are in average twice as large 
as those proposed bv lVaughan et al.l (12003). For each frequency, 
this procedure consists in randomly drawing a number of points 
N from the corresponding light curve, without avoiding repeti- 
tions and being the total number of data points in the original 
light curve. Then the Fvar parameter is calculated on the simu- 
lated light curve. This is repeated a statistically significant num- 
ber of times (100,000 here), in order to obtain a distribution of 
the Fvar values found. A Gaussian fit is then applied and the stan- 
dard deviation obtained is taken to represent the uncertainty of 
Fysi for the corresponding original data set. 

4.1 .2. Results on tine variability amplitude 

The fractional variability amplitude strongly depends on fre- 
quency, as shown in Fig. |4] Even though different in slopes, 

1 //eV 1 meV 1 eV 1 keV 1 MeV 1 GeV 




lO-'-'L , , , , , , , , , , , , , , , , J 

10' 10'2 10'= lO'S 102' 10" 
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Fig. 3. Average spectral energy distribution of 3C 273 from radio 
to gamma-rays (points), spanning from 4 to 44 years of observa- 
tions depending on the wavelength. The error bars are calculated 
as the standard deviation from the mean values and the grey area 
represents the observed range of variations. 
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Fig. 4. Fractional- variability, Fyai , spectrum of 3C 273 from ra- 
dio to gamma-rays (filled circles). Open triangles have been ob- 
tained after removal of flares in the infrared band. The dashed 
lines in the radio-millimeter, optical-UV and X-ray bands repre- 
sent the best fits. 



the radio-mm, optical-UV and X-ray Fyar all show a clear in- 
creasing trend. In particular the increase of Fyai is linear with 
the logarithm of the frequency with a slope of 0.287 + 0.004 and 
0.14 + 0.01 in the radio-mm and optical-UV bands, respectively 
(Fig.|4]i. In the X-rays, a single linear increase from 0.1 to 500 
keV does not represent the data well. There seem to be rather two 
different behaviours below and above ~ 20 keV, with a steeper 
slope of 0.11 + 0.02 at higher energies when compared to the 
much flatter trend in the lower band (0.03 + 0.01). 

The steep increase of Fy^ from radio to millimeter is well un- 
derstood in the frame of the evolution of synchrotron flares that 
are believed to pro pagate from highe r to lower energies, with in- 
creasing duration ( Tiirler et al.1200 0'). As a consequence, consec- 
utive flares overlap, creating an intense, continuous radio emis- 
sion with relatively small variations, whereas in the millimeter 
band very short flares stand out on a lower quiescent emission, 
resulting in a large fractional variability amplitude (Fig.|2]i. 

The position of the Fyar peak in the mm band, around 0.35 
mm (8.6 X 10" Hz), is uncertain due to the scarcity of data in 
the far-IR range. Nevertheless, a clear decreasing trend of Fy^ is 
visible in the IR band. In the 1-10 fim light curves (J to N) where 
the synchrotron flares have been subtracted, Fy^i- decreases by a 
factor of 2 or more compared to the same data sets including 
the flares. In particular, the variance of the M curve is smaller 
than the squared measurement uncertainties, whereas not enough 
points are left in the Q light curve to estimate Fyar- In the near- 
IR around 4 fim (L band), Fy^ again starts to increase towards 
shorter wavelengths, still staying below 0.1. The less variable 
component arising in the mid-IR band is likely to be thermal 
emission from dust. As it approaches the sublimation tempera- 
ture, the dust emission is contributing decreasingly, therefore an- 
other component has to be responsible for most of the observed 
radiation and variability at wavelengths around 1 yum (J and H 
bands). 

The fractional variability amplitude in the optical-UV do- 
main steadily increases from 0.1 to 0.2. This trend is consistent 
with the decomposition suggested by Paltani et al. (1998b), in 
which the red ('??') component, dominant in the optical range, 
shows smaller variations than the blue ('S') component, which 
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Fig. 5. Average mm to UV spectral energy distribution. The best 
fit model to the mm-IR data (continuous line) is obtained adding 
a cut-off power law (dashed line), a grey-body dust contribu- 
tion (dotted curve) an d the broken-power-law H component from 
iPaltani et all (Il998bh (dotted lines). In the optical-UV band the 
continuous line is the sum of the 'R and S (dot-dashed line) com- 
ponents. The mm to optical fluxes refer to the light curves where 
synchrotron flares have been subtracted. 



is dominant at UV frequencies, the total emission resulting from 
changes in the relative fluxes of these two components. 

Considering the continuous trend of Fyar from 8 x 10'^ to 
3 X 10'^ Hz, !R could extend towards long wavelengths to the 
near-infrared domain and be responsible for most of the emis- 
sion and variability properties in this range, outside the short- 
flare periods. The synchrotron emission, as extrapolated from 
the mm band, and the dust explain the far- and mid-IR emis- 
sion, but a cut-off or a break is required in the synchrotron ra- 
diation to avoid overestimating the near-IR flux. We thus fit the 
mm to IR spectrum (excluding synchrotron flares) with an expo- 
nentially cut-off power law (k^g^ y-Q ' mmg-- -) ^jj^j fj^ the contri- 
bution as modelled by Palta ni et al.l ( Il998b ). Paltani et al. found 
that the optical-UV data were consistent with H having either a 
single (= '??max case, with a slope of -0.7) or a broken power 
law shape (= 'R^in case, with slopes of 0.0 and -2.9). We choose 
to model 'R here with a broken power law with slopes of -0.4 
and -1.5, before and after the break, as this represents an inter- 
mediate shape between the ??min and ??rnax extreme cases. The 
fit results in a reduced chi-sqnaie x'^^^ - 1-5, mainly due to the 
deviation of the near-IR points from the model. Adding a com- 
ponent to account for the dust emission improves significantly 
the fit iY^^^ - 0.36, F-test probability - 2 x 10""*) and gives 

a„„ = 0.74+0.06, fcmm = 1.9x10-" at 10" Hz and Vc ^ 5x10'^ 
Hz (Fig. |5]l. The dust contribution is represented by an isother- 
mal grey -bod}Q with fitted temperature Tdust - 1156 + 92K and 
emitting area Adust = 1 -5 ± 0.5 pc^. The other grey-body parame- 



' An isothermal grey-body at temperature T emits at the observed 
frequency Vobs a flux density given by: 

f v'"(>'obs) = (1 + z) Ad„,, - e-^") Br(ve„, T) (2) 

where fi™(Vem, T) is the Planck function for a black body at temperature 
T, Di is the luminosity distance, A^^^t is the area of the projected source, 
Ty is the optical depth that can be expressed as a function of the dust 
emission index p and the frequency at which the optical depth equals 
1, Ty = (v eralvT=\f ■ TMs model does not take in to account the source 
geometry jPolletta et alfcOOdlTurler et al.ll2006l) . 



ters (i.e. the dust emissivity index and th e frequency at whic h the 
optical depth equals 1) are chosen as in iTiirler et al.l (12006 !). We 
apply the same model to the spectrum that includes synchrotron 
flares and to that derived using only synchrotron flares and we 
find that the same model with very similar parameters can rep- 
resent the different states of the mm to IR emission. The small 
differences between the T auhx an d Adust values we find and those 
reported by Tiirler et al] (|2006|) could be due to variability of 
the dust component, but can also be explained by the strong de- 
pendence of the dust modelling on the extrapolation of the mm 
emission which is quite uncertain due to the large variability of 
this energy band. 

In the X-ray domain, there is an increase of the fractional 
variability amplitude from 0.2 to 500 keV. However, the 0. 1-0.5 
keV curves are extremely sensitive to the addition or subtraction 
of a few data points, making highly uncertain the Fvar trend in 
this band. We can exclude the variations below 10 keV being due 
to the hydrogen column density, A^h, because this would imply 
extremely larger variations at low energies than are observed. 

For some gamma-ray light curves Fyai- cannot be calculated 
because of the small amount of available data or the large mea- 
surement uncertainties (1-10 MeV, 3-10 GeV). The increase of 
Fvar observed in the X-rays continues up to the gamma-rays, pos- 
sibly even above 100 MeV, because the observed decrease in the 
GeV range is most likely due to very large flux uncertainties. In 
addition, caution should be used when considering the variabil- 
ity results in the gamma-ray band above 30 MeV, because these 
light curves include also fluxes obtained from the same observa- 
tions by different analyses (see Sect. 12.5b . 

4.2. Characteristic time scales of variability 
4.2.1. Method 

A method to estimate the characteristic time scale of a light curve 
is given by the structure function (SF). For a tim e series F{t), the 
SF is defined as a function of the time lag t by jSimonetti et al.l 
Il985h : 

SF(t) ^< (F(t + t) - F(t)f > . (3) 

The SF is initially flat at a level given by the average measure- 
ment uncertainty, 2 ■ e^^, it then increases up to a maximum 
variability time scale, Tmax , indicating that the source variations 
are correlated. For r > T^ax, the SF flattens again at a value cor- 
responding to 2 ■ (S^ + e^^), where is the variance of the data 
(see an example in Fig.|6]l. 

We apply the structure function analysis to all data sets with 
at least 200 points and then further select those for which the 
upper plateau is clearly present, allowing a measurement of Tmax 
(Fig. IT). In a few cases, for the 0.8 mm and some of the optical 
light curves, the SF keeps increasing beyond the 2 • (S^ + Sen-^) 
value, indicating a r^ax longer than the sampled time scales (i.e. 
> 9 years). We exclude the light curves that have SFs for which 
it is difficult to reliably estimate the position of the plateau and 
the corresponding Tn,ax- 

In order to estimate the uncertainty on Tmax, for each data 
set we simulate 1,000 light curves (with a bootstrap procedure), 
compute the SFs and calculate Tmax- To automatize the process, 
we define Tmax as the shortest time scale for which the SF is 
higher than 97% of the expected value of the upper plateau. Then 
the distribution of the simulated Tmax values is calculated and 
the uncertainty on the original Tmax is calculated as the, possibly 
asymmetric, Icr standard deviation from the median value. 
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Fig. 6. Example of structure function for the 8 GHz light curve, 
renormalised to the expected value of the upper plateau, (2 • (S^ + 
e^^), continuous line). The dashed line is the expected value of 
the lower plateau (2 -e^^). The maximum time scale of variabil- 
ity is found to be around 4.4 years (vertical dotted line), when 
the SF flattens reaching the upper plateau. 



t 



I- 



109 



ia'2 



10'5 

Frequency [Hz] 



Fig. 7. Spectrum of the maximum variability time scale r^ax 
estimated from the structure function analysis (the lower lim- 
its are not displayed). The error bar calculation is explained in 
Sect. 14.2.11 The dashed line represents the best fit to the radio- 
mm data using Eq. (|4|. 



4.2.2. Results on the variability time scales 

The maximum time scales of variability from radio to mm fre- 
quencies are interpreted as due to the typical evolution of syn- 
chrotron flares, showing longer duration. The 22 and 37 GHz 
values of r^ax, calculated with a structure function analysis, have 
been already reported by iHovatta et al.l ([2007). Assuming that 
the measured variability time scales in the radio-mm band rep- 
resent an upper limit to the cooling time of electrons through 
synchrotron emission: 



Tmax(v) > fcool ^6x10*^ B 



D 



-1/2 



(4) 



(where fcooi is in seconds, B in Gauss, v in MHz and D is the jet 
Doppler factor), we can obtain a lower limit to the magnetic field 
responsible for this emission. This is just a rough lower limit of 
B as Eq. (|4]i is obtained in the case of a plasma cloud cooling 
through synchrotron mechanism in a constant magnetic field, a 
rather different condition than that of a shock in a jet. The fit of 
Log(Tinax) vs Log(v) with Eq. (gll results in B > Bmin - 0.089 + 
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Fig. 8. X-ray structure functions in the 9-20 and 20^30 keV 
range, based on RXTE (circles) and BATSE (squares) data, re- 
spectively. Both SFs have been renormalised to the expected 
value of the respective upper plateau (continuous line). The max- 
imum time scales of variability are found to be around 0.22 and 
1.3 years for the lower and higher energies, respectively (vertical 
dotted lines). 



0.003 Gauss (Fig. ID, when D = 5 dSavolainen et al.ll2006l and 
references therein). This result is consistent with the magnetic 
field range found by previous studies, where an estimate of B ~ 
0.4 - 1 Gauss was obtained using the time scales of sub-mm to 
cm and IR flares (Courvoisier et al. 1988; Robson et al. 1993). 

In the UV range, Tmax is basically con stant (as already found 
bv lPaltani et al.lll998bl;lFavre et al.ll2005h and much shorter than 
in the optical range, where Tmax is about 4 years or even longer 
than the sampled time scale s. This is due to th e dominance in the 
UV of the S component (Palt ani et aT]|1998bl) . 

In the X-rays, a significant difference can be seen in the max- 
imum variability time scales below and above 20 keV. In fact, the 
BATSE 20-70 and 70^30 keV SFs reach the upper plateau at 
a time r^ax = 1 -3 + 0.4 years, hereas in the 4-9 and 9-20 keV 
bands the flattening of the SFs happens at r^ax - 0.219 + 0.003 
years. The difference in the measured maximum time scales can 
be clearly seen in Fig. [8] where the 9-20 keV and the 20^30 
keV (from combining the two BATSE data sets) SFs are com- 
pared, on time scales from 1 day to 10 years. Even though the 
longer-timescale variations of the BATSE light curve occured 
before 1996, i.e. outside the period covered simultaneously with 
RXTE, the 50 and 100 keV SFs (not containing BATSE data, but 
including RXTE data) are consistent with BATSE results, show- 
ing a smooth flattening in the 0.5-2 years range. 

4.3. Cross-correlation analysis 

We correlate all the data sets overlapping for at least 4 years and 
with a relatively dense sampling during this period, concentrat- 
ing on the correlations with the X-rays where this analysis is 
possible for the first time. We do not discuss here optical-UV 
correlations, already presented bv lPaltani et al.l (Il998bl) and for 
which we obtain similar results. 

Two methods can be used to estimate the correlation func- 
tion: the interpolated correlation analysis (ICF) where both sig- 
nals are independently interpolated, and the discrete correlation 
analysis (DCF), where the correlation function is calculated in 
several time bins. Both methods have strengths and weaknesses 
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Fig. 9. ICF cross-correlation functions showing the delay of the 
37 GHz light curve on some representative radio-nun data sets at 
5 GHz, 8 GHz, 15 GHz, 37 GHz (auto-correlation) and 0.8 mm. 

depen ding on the data and their sampling (IWhite & PetersonI 
11994 . We apply both the ICF and the DCF methods (u sing 
the improvements prop osed by iGaskell & PetersonI 1 19871 and 
I White & PetersonI 1 1994i) and use the consistency between ICF 
and DCF as a criterion for the robustness of the correlations. In 
addition, we estimate the uncertainties on the correlation func- 
tions using a bootstrap procedure. For better clarity, only for one 
DCF (top left panel in Fig. fTTT i and one ICF (centre left panel 
in Fig.fTTTi we show the error bars in the Figures. The time lags 
are corrected for the effect of redshift. We discuss separately in 
the following the correlation of the medium (< 20 keV, mainly 
RXTE data) and hard (20-70 keV, BATSE data) X-ray curves 
with the other wavelengths. 

4.3.1. Significant correlations 

- Radio and mm: the radio and mm light curves are very 
well correlated, showing peaks with correlation coefficients 
R > 0.6 + 0.01. The observed delays indicate that the higher- 
frequency emission leads the lower-frequency one (Fig. |9]), 
in agreement with what is expected from the evolution of 
sync hrotron fla res that propagate from higher to lower ener- 
gies jTtirler et al. 1999a, 2000). 

- IR and optical- UV: when the synchrotron flares are excluded 
from the data sets, a main peak is observed in the correlation 
between these light curves (Fig.fTOt. with correlation coeffi- 
cients reaching (0.5-0.8) +0.04 at lags between -1.2 and -1 
years for the K curve and shifting to lower lags for H and J 
(between -1 and -0.8 years for the latter). This indicates that 
the IR emission follows the optical-UV one with a delay that 
increases with IR wavelength. There is no clear dependence 
on the optical-UV frequency of the time delays and shape of 
the correlation function, apart for a possible small-lag peak 
(at -0.3 to years) for the J curve which is more prominent in 
the optical than in the UV correlations. It is therefore proba- 
ble that the IR correlates with both the and S components 
rather than with only one of them. 

- Hard X-rays and radio-mm: there is a correlation between 
the 5-37 GHz light curves and the hard X-rays, with the 
higher radio frequencies leading and the X-rays correlating 
at about lag with the 5 GHz emission (centre left panel in 
Fig. ED. The correlation peaks reach R = (0.5 - 0.6) + 0.06 
and are quite broad, so that a precise time lag is difficult 
to estimate. The cross-correlation plots with the mm band 
show a significant peak at ~ -2.5 years (centre right panel 
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Fig. 10. Cross-correlation functions (calculated with the ICF 
method) for the near-IR light curves K, J and H with some of 
the optical and UV data sets. In a correlation 'X'-'Y\ a posi- 
tive time lag means that Y is leading. Top panel: V - IR. Middle 
panel: 3000 A - IR. Bottom panel: 1300 A - IR. 



in Fig. [TTi) . which would indicate that the X-rays are lead- 
ing. The radio-mm correlation with the hard X-rays is due 
to the matching of the single, strong peak characterising the 
BATSE light curve with one synchrotron flare (Fig. |2]), i.e. 
the flare propagating in 1991-1993 from 37 to 5 GHz and the 
1997 flare in the mm. This strongly questions the credibility 
of the correlations, and therefore, a more complete theoreti- 
cal understanding is needed to know whether this correlation 
is due to chance occurrences or whether there is a connection 
between the hard X-ray and the jet emission and a common 
mechanism is responsible for both emissions. 
- Hard X-rays and optical-UV: there are two peaks in the cor- 
relations between these bands (see lPaltani et al.l ll998a'). one 
at ~0 lag and the other at about 1 .5 years, with the blue bump 
emission leading the X-rays. The peak without delay is dom- 
inant in the optical-X-ray correlation (bottom left panel in 
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Fig. 11. Cross-correlation functions as calculated with the ICF 
(continuous lines) and the DCF (circles) methods. In a correla- 
tion 'X'-'F, a positive time lag means that Y is leading. Top 
left panel: 20-70 keV - 5 keV, with error bars on the DCF. Top 
right panel: 5 keV - 15 GHz. Centre left panel: 20-70 keV - 5 
GHz, with error bars on the ICF. Centre right panel: 20-70 keV 
-1.3 mm. Bottom left panel: 20-70 keV - V. Bottom right panel: 
20-70 keV - 3000 A. 

Fig.fTTTl and the peak at 1.5 years is dominant in the UV-X- 
ray one (bottom right panel in Fig.fTTI). This seems to indi- 
cate that the H component, dominant in the optical but still 
contributing in the UV, correlates at lag with the X-rays, 
whereas the S component is emitted 1.5 years earlier These 
correlations, as those between the radio-mm emission and 
the hard X-rays, are driven by the strong peak of the BATSE 
light curve. However, the less structured light curves in the 
optical (compared to the flaring structure of the radio-mm 
ones) and the fact that we find exactly a 0-lag correlation 
make the optical-hard X-ray relation more reliable. Since 
the UV-hard X-ray correlation shows a similar shape to the 
optical-hard X-ray one, the reliability of the former is likely 
related to that of the latter. 



4.3.2. No correlations 

- Medium and hard X-rays: this correlation function shows 
a few peaks and repeated oscillations (top left panel in 
Fig.fTTI). Two main peaks are obtained with the ICF method 
atOand-1 yearwithT? ~ (0.5-0.6)+0.1, whereas in the DCF 



12 



a. 

" 10 



20 22 24 26 28 30 32 
3000 A flux [mjy] 

Fig. 12. 2 keV versus 3000 A flux obtained from quasi- 
simultaneous (within one day) observations. The Spearman-rank 
test results in a correlation coefficient of R - 0.25 and a proba- 
bility for chance occurrence of f = 0.3, indicating that there is 
no correlation between these two data sets. 



these peaks are at ~ (0.3 - 0.4) + 0.15 and an additional 
peak appears at lags of -2.5 years, with R ~ 0.45 + 0.2. The 
discrepancies between the ICF and DCF and the large error 
bars of the DCF imply that any correlation is questionable. 
Medium X-rays and radio-mm: no correlation is found (top 
right panel in Fig.fTTTl. suggesting that the main medium-X- 
ray variability is not related to the jet flaring emission. A cor- 
relation between the IR and the medium X-rays during two 
flares has been interpreted as evidence for an SSC origin of 
this X-ray emission (McHard v et al.ll2b07h . The lack of cor- 
relation between the medium X-rays and the radio-mm sug- 
gests t hat the SSC mechanism identified bv iMcHardv et al.l 
(|2007|) would rather be a secondary effect. 
Medium X-rays and optical-UV: as the overlap between 
these data sets is not sufficient for a correlation analysis on 
time scales of years, we concentrate on quasi-simultaneous 
(within one day) optical-UV and X-ray observations. We ap- 
ply a Spearman-rank test to quantify the correlation in these 
subsets of data, containing between 6 and 50 data points. 
Apart for one significant correlation between 1 keV and 3000 
A (correlation coefficient R = 0.7 with a probability for 
chance occurrence P = 0.002) and a marginal one between 
5 keV -B (R - 0.4, P - 0.03), no other correlation is found 
(see Fig. [T2]for an example). This result questions the con- 
nection between optical-UV and X-ray emission which is ex- 
pected in reprocessing models (see Sect. 15. 3t . 



4.4. X-ray spectral evolution 

We study the evolution of the X-ray photon index a obtained 
from spectra extracted within the 2-20 keV range, to exclude 
the soft-excess variations and the hard-energy domain. We con- 
centrate our attention on the last 10-years d ata (1996-2006) as 
they r epresent the best-sampled period (see IChernyakova et all 
f2007f for a discussion of a variations over 30 years). 

The 2-20 keV photon index has an average value of 1.71 + 
0.07 and exhibits significant variations of about Aa - 0.2 on 
short time scales (Fig.fTSJ. It reached a peak of apeak — 

1.84 dur- 
ing 2003 and since then the spectrum has been hardening again. 
This feature dominates the evolution of the spectral shape of the 
last 10 years and it does not seem to correspond to any of the 
trends observed in the light curves of our database. In fact, no 
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Fig. 13. Evolution of the X-ray, 2-20 keV, photon index between 
1996 and 2006. RXTE data are rebinned to one week for better 
clarity. 



correlation is found with the light curves that significantly over- 
lap the 1996-2006 photon index curve, i.e. the radio, 3.3 and 
0.8 mm and the X-ray ones. We investigate in more detail the 
relation between the 2-20 keV photon index and the 2-20 keV 
flux. Depending on the time range considered, there is no corre- 
lation (e.g. 1996-2006), correlation (1999-2000, Spearman cor- 
relation coefficient R - 0.55, P = 3 x 10"^) and anti-correlation 
(2001-2005, R = -0.26, P = lO""*) between X-ray photon in- 
dex and flux. These relations appear to represent random oc- 
currences rather than a real photon index-flux connection. If the 
1999-2000 correlation were the signat ure of the accr e tion d isc 
or corona component, as suggested by iKataoka et al.l (|2002|) . a 
similar trend should be observed in 2001-2005 data when the 
jet emission reached its historical minimum (Tiirler et al. 2006) 
and the disc/corona component should have dominated the X-ray 
spectrum. Instead, an anti-correlation is observed in 2001-2005, 
in contradiction with this interpretation. 



5. Discussion 

5. 1 . The synchrotron flaring emission 



iTiirler et"aLl (1 1 999ai l2000h have extensively studied the variabil- 
ity behaviour of 3C 273 in the radio-mm domain and shown 
that it can be well reproduced with the superposition of syn- 
chrotron flares from shock waves, starting as short flares at high 
frequencies and propagating to lower frequencies while evolv- 
ing on longer time scales. This scenario explains the decreasing 
amplitude of the variations (Fig. [Hi, the increasing characteristic 
time scales (Fig.|7j and the correlations (Fig.|9]l from the mm to 
the radio band. It is likely that some flares start as very short and 
peake d spikes in the IR or even visible bands ( iCourvoisier et al.l 
[l988h . 

5.2. The origin of the "R component 

The two components "R and S proposed bv lPaltani et al.l (Il998bh 
to model the blue bump show different spectral and variabihty 
properties, suggesting different origins of the underlying emis- 
sion processes. 



In view of the measurement of slightly polarized emission 
in the optical (even outside the periods of synchrotron flares, 
|de Diego et al. 1992), one natural explanation is that R could 
be synchrotron emission not directly related to the radio-mm jet 
(Fig.|5]l. The small amplitude of the variations observed in the H 
component (Fig.|4]l could be explained by an electron population 
with cooling time longer than the injection time in the emitting 
region. Even th ough is no t compl etely constrained by the de- 
composition of lPaltani et al.l (Il998bl) . its global spectral shape is 
also compatible with optically thin synchrotron emission. 

We investigate also the possibility that H origins in the broad 
line region (BLR). Considering that optical emission lines are 
estimated to contribu te at most by 10% of the total continuum 
(iPaltani et al.lll998bl) and an additional 10% is expected from 
the B akner continuum and the Feii emission (' Courvoisier et al.l 
119871) . the remaining "R emission could have a thermal origin, as 
Bremsstrahlung radiation from the BLR. This possibility is sup- 
ported by the fact that, even though most of the flux variations of 
Lya and Civ has been found to correlate with the photoionising 
UV emission, the long-term variability properties of these lines 
are m ore similar to those of the R com ponent (jP altani & Tiirle^ 
l2003h . From photo-ionisation models dRobson [l996l) . we cal- 
culate a BLR mass o f 160 Mq, using an average Ha flux of 
6 X 10" erg cm s-i dKaspi et al.ll2000l) and assuming uniform 
electron density of the BLR with a canonical value of 10'" cm"^. 
If the BLR gas is clumpy, the BLR mass is expected to be less 
than an order of magnitude larger. A mass of 160 Mq at 3000 K, 
as suggested by the spectral shape, results in a Bremsstrahlung 
luminosity of LBremss = 1 .7 X lO"'^ erg/s, which is only 1 % of the 
total luminosity observed in the R component (4 x 10"*^ erg/s), 
ruling out a BLR Bremsstrahlung origin for R. 



5.3. S and X-rays relation 

A short-lag correlation between the UV and X-ray (< 10 keV) 
emissions is expected if reprocessing is the main process pro- 
ducing the UV radiation. In the p ast, such a correlation has been 
found dChernvakova et alj l2007b or not ( Walter & CourvoisTeJ 
19921) on smaller data sets, depending on the time period. We 
obtain here that there is no short-lag correlation between the 
light curves on these two bands, indicating that, if reprocess- 
ing is the driving proce ss (as sugge sted by the UV spectral and 
variability properties, Paltan i et al.|[T99 8b). somehow the X-ray 
radiation interacting with the UV source is not the same as that 
observed from Earth. This could happen if the X-ray emission 
IS anisotropic (IUttlevll2005h or if the X-ray source heating the 
disc is positioned very close to the black hole. In this case, only 
a small part of the X-rays could escape, whereas most of them 
would be gravitationally bent by the black hole and illuminate 
the accre tion disc. This possibility based on the light bending 
model of iMiniutti & FabianI (12004 ') needs proper simulations in 
order to prove that this geometry would actually reproduce the 
observed emission. 

Several studies on Seyfert galaxies have given conflicting 
results on th e optical/UV-X-ra y relation, showing correlation 
(NGC 5548, Uttley et al."2003') or no correlation (NGC 3516, 
iMaoz et al. 2002), sometimes depending on the observed time 
scale (NGC 4051, Shemmer et al.l2003b . and with different vari- 
ability amplitude in the optical/UV compared to the X-rays (see 
lUttlevI [20051 for a review). There are however several reasons 
for which physically connected blue bump and X-ray emissions 
would not be correlated. For example if the X-rays are due to 
Comptonisation of the UV photons, the UV-X-ray flux relation 
could depend on a third parameter, e.g. the X-ray spectral index 
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(lNandrall2001l ). Yet, e ven though such a relation wa s detected in 
3C 273 in the past bv lWalter & Courvoisieii ( IT992I) . when more 
recent data are added no correla tion of these three parameters is 
found (IChernvakova et af]|2007h . Among other possibilities, the 
feedback between reprocessing and Comptonisation or the ge- 
ometry and the size of the X-ray source or of the reprocessing 
medium could wash out the variations of the observed emission 
dlandra 2001). 

An alternative scenario for the production of UV radia- 
tion involves n iatter accreted in clurnps, ra ther than through an 
accretion disc dCourvoisier & Tiirleil 12005 !). The interaction of 
these clumps at ~ 100 Schwarzschild radii generates optically 
thick shocks producing the UV emission, whereas optically thin 
shocks closer to the black hole give origin to the X-rays. In this 
model, both a short-lag and a longer-lag correlation are expected. 
The latter correlation has been found in 3C 273 but between the 
hard X-rays and the UV. The application of this model to the 
case of 3C 273 will be the subject of a forthcoming paper. 



5.4. The near-IR emission 

The amplitude of the variations (Fig. |4|i and the shape of the 
near-lR to UV spectrum (Fig. |5]l suggest that fi could extend 
to near-lR wavelengths and be responsible for the variability 
in this band, synchrotron flares excepted. An additional dust 
contribution has b een proposed in the past (iRobson et alJll983t 
lTurleretalJl2006l) and is confirmed here by the need of a ther- 
mal component to model the average near-IR spectrum. As the 
extrapolation of the mm emission is higher than the near-IR flux, 
a cut-off' is required in the high frequency synchrotron emission 
at ~ 5 X 10'^ Hz, and therefore H and dust emission are likely 
the dominant contributors in the near-IR. 

The interplay between the components is observable in the 
correlations between the near-IR and the optical-UV emission 
(Fig. [Tot . In the most likely scenario, the dust is heated up by 
the UV/S continuum. The complexity of the dust component 
appears in the 0.8-1.2 y delays found in the correlations. These 
could represent the travelling time of the S photons from the 
UV source to the dust location. The delays increase with in- 
creasing IR wavelength, which is expected from an extended 
dust region with a distribution of temperatures. A similar cor- 
relation between near-IR and optical/UV has been observed in a 
number of Seyfert 1 galaxies. The time lag of the K with the V 
band seems to be proportional to the square root of the optical 
luminosity dSuganuma et al.l2006l and references therein) and is 
believed to represent an estimate of the inner radius of the dust 
torus in these objects. Even though 3C 273 is 100 times more 
luminous than the Seyfert 1 Fairall 9, their time lags are similar 
dClavel et al.lll989h . indicating that 3C 273 does not foUow the 
lag-luminosity relation of Seyfert 1 galaxies. 

The correlation of the J light curve with the optical-UV 
shows a secondary peak at short lags, which is higher in the opti- 
cal than in the UV band (Fig.fTOli. The same effect, even though 
much weaker, is still present in the correlation of H with the 
optical-UV. This secondary peak could be the signature of the 
K component, dominant in the optical and responsible for most 
of the observed emission also in the J band. The same effect 
has been reported for Fairall 9, for which the J band was found 
to be correlated without lag to the optical-UV curves, whereas 
K a nd L were correlated with delays of ^ 1 year (Clavel et al. 
11989 ). This correlation represents an additional evidence that fi 
contributes to the near-IR emission of 3C 273. 



5.5. The origin of tlie X-ray emission 

In the X-rays, the maximum variability time scales are signif- 
icantly longer above 20 keV than below and the amplitude of 
the variations is increasing steeply with energy above 20 keV, 
whereas at lower X-ray energies the increase is rather shallow. 
In addition, no convincing correlation is detected between the 
light curves below and above 20 keV. 

The difference between medium and hard X-rays imply the 
presence of two varying parameters, which can be interpreted 
either as two distinct physical components or as a single compo- 
nent with varying spectral shape. Both descriptions are equiva- 
lent, they may, however, lead to different scenarios. 

In the two-components scenario, the high-energy emission of 
3C 273 could be due to the superposition of a Seyfert-like and a 
blazar-like component, as suggested by e.g. iGrandi & Palumbol 
(2004). Comptonisation of thermal plasma for example in an ac- 
cretion disc would produc e X-ray emissio n similar to that ob- 
served in Seyfert galaxies dGrandi & Palum bo 2004) . The occa- 
siona l detection of a iron Kg emission fine (e.g. iKataoka et al] 
I2002t IGrandi & Palumbol l2004t IChernvakova et aljr2007h con- 
firms the similarity between 3C 273 and Seyfert emission in 
the medium X-rays. Above 10-20 keV, inverse Compton on a 
population of relativistic electrons from a non-thermal plasma 
would produce the observed high-energy emission, as suggested 
by the presence of a break in the sp ectrum around 1 MeV 
dJohnson et al.lll995HLichti et al.ll 19951) . 

In the single-component scenario, the 2 to 500 keV spec- 
trum could be produced, for example, by a single power law 
pivoting at low X-ray energies. Then the variations at medium 
X-rays would be dominated by the variation in the low energy 
flux of this component, whereas the variations at hard X-rays 
would rather be due to the variation of the spectral index. The 
long time-scale variations of the spectral index (Fig.[T3]l and of 
the hard X-ray emission (Fig.jSJ could support this decomposi- 
tion. 

The lack of short-lag correlation between the mm and the 
X-ray radiation indicates that there is no clear connection be- 
tween the X-rays and the radio-mm flaring emission from the jet. 
Compton losses are thought to be the dominant electron cooling 
process in the first phase of the evolution of the shocks in the jet 
(iMarscher & Geadl985h . Therefore, if this process were produc- 
ing the bulk of the X-ray emission, the mm flares would follow 
the X-ray ones with short delays, which is not observed. In addi- 
tion, the hard X-ray peak in the 20-70 keV light curve is broad 
(Fig. |2]i, also not expected if it represented the early phase of 
the shock flaring emission. The lack of correlation between the 
medium X-ray and the mm suggests that the simultaneous IR 
and medium X-ray emission detected by McHardy et al. (20O3) 
and interpreted as due to SSC from the jet must therefore be the 
result of a secondary process. This shows that some synchrotron 
events can be so energetic that they dominate the X-ray emission 
even up to 20 keV for short time periods. 

A correlation without delay is observed between the hard 
X-rays and the 'R component. If H is synchrotron emission, as 
suggested here, the hard X-rays could be produced by the same 
electron population emitting the 'R/synchrotron component. An 
alternative to this synchrotron self-Compton mechanism (SSC) 
is an external Compton process (EC) suggested by the 1 .5 y de- 
lay observed between the S emission and the hard X-rays. If one 
assumes that S provides the seed photons for the Compton pro- 
cesses producing hard X-rays, this 1.5-year delay could indicate 
the distance covered by the photons from the UV to the X-ray 
source. 
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6. Summary and conclusions 

We have presented an update of the multiwavelength 3C 273 
database with the last 10 years of data and performed an analysis 
of the variability properties of this object with different methods. 

In order to explain the significantly different timing proper- 
ties of the X-ray emission below and above ~ 20 keV, either two 
separate components (possibly a Seyfert-like and a blazar-like 
one) or at least two parameters with distinct variabiUty proper- 
ties are needed. The 5-500 keV emission does not seem to be 
related to the mm-jet, questioning the usually assumed common 
origin for the electrons emitting the radio-mm and the hard X- 
ray radiations. On the other hand, the electron population at the 
origin of the dominant hard X-ray emission could be the same as 
that responsible for the optical emission, which is likely to orig- 
inate from synchrotron mechanism. If the Comptonised photons 
are the U V ones, the 1 .5-year delay between the U V and the hard 
X-ray emission could be an indication of the distance between 
the X-ray emitting region and the seed photon source. 

The blue-bump emission does not seem to be related in an 
obvious way to the X-ray emission below 20 keV. If there is a 
close interplay between the blue bump and the X-rays, as ex- 
pected in Comptonisation and reprocessing models, then some- 
how the medium-X-ray emission we observe is not the same as 
that interacting with the blue bump, maybe due to some partic- 
ular geometry of the emitting region or anisotropy of the X-ray 
radiation, possibly because of Ught bending effects. 

The near-IR spectral and timing properties are well ex- 
plained by the combined contribution of emission from dust il- 
luminated by a UV source at < 1 light year and the extension of 
the optical/synchrotron component. 

The variabiUty properties of the radio-mm emission are well 
understood in the frame of synchrotron flares from shock waves 
propagating along the jet, with spectral and timing characteris- 
tics evolving with time. Assuming that the maximum variabiUty 
time scales measured in this band represent an upper limit to the 
electron cooling time via synchrotron emission, we obtain an es- 
timate of the magnetic field of B > 0.09 Gauss. 

The complexity of 3C 273 emission (and of AGN in general) 
is, at least in part, due to the presence of many emission compo- 
nents of different nature but similar luminosities, probably origi- 
nating in different locations within the AGN. In addition, in sev- 
eral energy domains the observed radiation is produced by the 
superposition of two or more of these different components. This 
highlights the importance of long-term multiwavelength studies 
for a better understanding of the physics of these objects. An 
easy access to all observational data is therefore crucial to be 
able to use them even a long time after the observations were 
performed. 
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Table 1. Radio to y-rays 3C 273 light curves used for the variabihty analysis. Col.l: light curve name. Epoch: the epoch between the first 
and last observation. A^: number of data points used for the variability analyses, v: average frequency. Fy and cr^: average flux and its standard 
deviation in Jy. vF^: average intensity in units of 10"'' erg cm"^ s~' . F^^: variability amplitude and its uncertainty from bootstrap simulations. 
Tmax: maximum variability time scale from structure function analysis and its uncertainty from bootstrap simulations. '. . .' indicates those 
light curves not used for the particular analysis due to different reasons, for instance an insufficient number of data points, (wf): data sets from 
which the synchrotron flares have been excluded. flux in counts/s. flux in 10""* erg cm"^ s"' . 
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'''max 
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■ ■ ■ 


10 GHz 


1965-1994 
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1983-1995 


17 
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6.6 


0.73 ±0.10 




Q 
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17 
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0.31 
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0.35 + 0.07 




N 


1970-2004 


75 
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0.35 
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0.23 ± 0.02 
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1970-2004 


54 
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0.31 
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0.14 + 0.04 




M 
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33 
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-2 
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-3 
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0.150 + 0.008 
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-2 
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-3 
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-2 
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-3 
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Table 1. continued. 



S. Soldi et al.: The multiwavelength variability of 3C 273 





Epoch 


N 
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Hz 


' V 
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0.1 keV 
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-1995 


20 
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2.7 10" 


-4 


7.1 10" 


-5 
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0.24 ± 0.06 


0.2 keV 
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17 


4.84 10^^ 


1.1 10" 


-4 


1.7 10" 


-5 
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0.16 ±0.03 


0.5 keV 


1979- 
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58 


1.21 10" 
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-J 
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-0 


3.8 


0.20 ± 0.02 


1 keV 


1978- 


-2005 


69 


2.42 10" 
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_5 
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3.9 


0.22 ± 0.02 


2keV 
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-2005 


93 


4.84 10" 


8.7 10 


_5 


2.6 10 




4.2 


0.27 ± 0.02 


5 keV 


1970- 


-2005 


1032 


1.21 10^^ 


4.7 10" 


-6 


1.2 10" 


-6 


5.7 


0.249 ± 0.006 


4-9 keV 


1996- 


-2005 


2567 


1.57 lO^** 


6.2" 




1.4" 






0.230 ± 0.003 


lOkeV 
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-2005 


1026 
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-6 


8.2 10 


-7 


7.0 


0.265 ± 0.007 


9-20 keV 


1996- 


-2005 


2567 


3.51 10^^ 


2.3" 




0.6" 






0.239 ± 0.004 


20keV 


1976- 


-2005 


987 


4.8410^^ 


1.810" 


-6 


5.4 10" 


-7 


8.7 


0.275 ± 0.008 


20-70 keV 


1991- 


-2000 


91 


1.09 10^^ 


2.6'' 




1.2* 






0.41 ± 0.04 


50keV 


1977- 


-2005 


1003 


1.21 10^^ 


9.6 10" 


-7 


3.610" 


-7 


11.6 


0.33 ± 0.02 


100 keV 


1978- 


-2005 


1001 


2.4210''' 


5.910 


-7 


2.3 10" 


-7 


14.3 


0.34 ± 0.01 


200 keV 


1978- 


-2005 


41 


4.84 lO^'' 


4.8 10 


-7 


3.2 10 


-7 


23.2 


0.40 ± 0.08 


70-430 keV 


1991- 


-2000 


91 


6.04 lOi' 


4.2* 




2.2* 






0.47 ± 0.05 


500 keV 


1990- 


-1999 


15 


1.21 102" 


2.8 10" 


-1 


1.710 


-7 


33.9 


0.43 ± 0.09 


IMeV 


1990- 


-1999 


20 


2.42 102" 


1.710 


-1 


9.9 10 


-8 


41.1 




30MeV 


1991- 


-1999 


17 


7.25 10^1 


3.3 10" 


-9 


2.5 10" 


-9 


23.9 


0.58 ±0.15 


100 MeV 


1976- 


-1997 


28 


2.42 10^2 


3.610- 


10 


2.610- 


10 


8.7 


0.64 ± 0.07 


300 MeV 


1976- 


-1997 


16 


7.25 10^2 


8.3 10- 


11 


5.5 10- 


11 


6.0 


0.48 + 0.13 


1 GeV 
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-1997 


15 


2.42 1023 


1.410- 


11 


9.4 10- 


12 


3.4 


0.20 ± 0.17 


3GeV 
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15 


7.25 1023 
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12 
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12 


2.1 
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1991- 


-1997 


14 
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13 
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13 
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yr 
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1.4 
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+0.003 
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+0.002 
-0.003 



+0.3 
-0.6 



